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CHRONOGRAPH  ERROR  IN  DRAG  MEASUREMENTS 

ABSTEIACT 

The  precision  of  the  least  squares  method  of  fitting  time -distance 
data  from  an  observed  trajectory  is  considered  for  the  case  where  time 
is  assumed  to  be  a  polynomial  in  distance.  A  relation  between  the 
variances  of  time  and  distance  errors,  assumed  in  a  previous  work,  is 
derived.  With  this  relation  and  a  model  of  digital  chronographs,  the 
variance  of  the  timing  error  is  found  for  selected  chronograph  arrange¬ 
ments;  in  a  favorable  arrangement,  this  veiriance  is  only  one-twelfth  as 
large  as  the  squared  counter  resolving  time.  A  comparison  of  time  and 
distance  errors  is  made. 
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I.  INTRODUCTION 


The  direct  free-fllght  technique  of  determining  the  drag  properties 
of  a  missile  requires  measurements  of  the  time-distance  history  of  a 
model  along  a  trajectory.  The  distance  measurements,  usually  obtained 
from  flash  photographs  of  the  missile  and  a  fiducial  marker,  are  com¬ 
bined  with  simultaneous  time  observations,  which  are  ordinarily  obtained 
from  digital  chronographs  coordinated  with  the  flash  pictures.  These 
observations  are  used  in  a  least  sqxoares  procedure  to  evaluate  the  drag 
properties  of  the  model.  The  accuracy  of  the  method  has  been  analyzed 
by  Karpov^,  with  special  attention  paid  to  the  spatial  arrangement  of 
the  observation  stations.  The  effect  of  the  pure  timing  error  inherent 
in  digital  type  chronographs  has  not  been  as  well  clarified;  and'from 
the  following  discussion,  the  size  of  this  error  appears  to  be  consider¬ 
ably  smaller  than  the  counter  resolving  time,  which  has  been  frequently 
used  as  a  measure  of  the  timing  error.  In  conventional  digital  chrono¬ 
graphs,  the  counter  resolving  time  is  necessarily  finite,  and  the  com'^ 
plexity  and  expense  of  the  chronograph  Increase  rapidly  as  the  resolving 
time  is  decreased.  An  appreciation  of  the  magnitude  of  timing  error  is 
necessary  for  fixing  the  specifications  of  the  chronograph  needed  for 
data  of  a  given  accuracy. 

II.  PRECISION  OF  THE  LEAST  SQUARES  METHOD 


II-l.  Polynomial  Fitting  ; 

It  is  assumed  that  the  time,  t,  when  the  missile  is  located  at  a 
particular  distance,  Z,  along  the  trajectory  is  related  to  the  distance 
by  an  algebraic  polynomial  of  degree  q,  so  tha,t 


q 

^  =  ZI 

J=0 


(1) 


The  a.  are  the  coefficients  of  the  polynomial  and  typify  the  motion. 
The  derivative  of  t  with  respect  to  z  represents  the  reciprocal  of  the 


missile  velocity,  V,  aloog-tfee  trajectory-,  thus-, 


i  =  vv=  z: 


j=o 


J  a.  z 

kj 


j-1 


(2) 
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From  tHe  free-fl'l^t  experiment,  discrete  pairs  of  values  of  time  and 
distance,  t^  and  are  ottained'at  Ndtservation  stations .  TEie  number 
N  is  usually 'arr^ged  "t6''be  Targe r  tlian  q  +  1,  the  number  of  uhSnown 
cdefficierits  ih'Eq(l),  in  order  that  a  least  squares'  techhique  Tl)e‘ ai^li- 
cable.  ' 

We’  define ’the  Tcth  residual,  R^,  to  be  the  distance  parallel  to 
the  t  axis ‘from  the  "''■bes't''  curve  "to  the  k'th  data  point;  or. 


^  ^  ^k 


(5) 


Fittin^'Toi’  the  smallest  value  of 'the  sum  of  the  squa'red'  residuals  in 
the  cohve'htibnaT  inanher  requires  that  the  a  -  satisfy"  the  noittaT  equations 

J 


q:  •<:,  CE  -  ■  V 


(4). 


Here  Z  is  a  q  +1  square,  symmetric  matrix,  with 


N 


Hij  =  ST 

k=l 


l+j 


(5) 


and  a  and  t  are  column  vectors  each  of  q  +  1  components,  as  foUowa: 


CoO,  =  a. 


and 


r  T  N  1 

Ml  -  C  4 


(6) 

(7) 


k=l 


'where  i  arid  j  range  "from "0  to  q.  A  speci'aT  case  where  the  distance 
data  is  symmetric  about  a  central  origin  is  riotahle  because  of  the 
vanishing  of  the  particular  elements  of  the  Z  matrix  where  i  +  j  is 
odd 
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II-2.  "Erfect  of  Time  Origin 


The' addition  of  an  arbitrary  constaJit  b^  to  each  t^^  affects  only 


the  value  of  a^  in  the  a  vector.  This  result  is  seen  by  observing  that 
the  addition  of '  a  vector  (whose’ Urst  component  is  b^  and  all  other 
components  are  zero)'  to  the'O!  vector  makes' the  product  Z  (o;  +  p)  equal 


to  T  +  Tq,  where  the  components  of  are  just  those  of  the  first 
column' of  ' the  'Z  matrix  multiplied  by  b^j  but  t  .+  also  resiilts  if  the 
value  “Bq 'is  added  to  each  t^^  in  the  t  vector.  'This  fact  is  not  unexpected 
(since  it  merely  corresponds  to  a  simple  change  of  time  origin  which  will 
not  affect  the  intrinsic  properties  of  the  polynomial)  and  proves  useful 
in  the  later  discussion. 

II-3-  Error  of  Fit 

To  investigate  errors  in  the  a.  due  to  errors  in  t  and  z  ,  it  is 

i  iC  iC 

observed  that^  although  the  z  values  may  be  in  error,  the  least  squares 

■  ■  iC 

fit  is  made  with  residuals  in  the’ t  direction  only.  However,  the  con¬ 
tribution  of  the  distance  errors  should  be  included  when  evaluating  the 


precision  measure  for  the  a  components.  Since  a  is  a  function  of  both 
the  and  the  t^,  we  have: 


N 

da  -  ^ 


Sa 


N 


(8) 


The  variance  of  da  is  found  from  the  sum  of  the  variances  of  the  two 
right  hand  members  of  Eq(8),  since  the  error  distribution  in  time  and 
distance  are  independent.  A  limiting  upper  bound  for  da  can  be  found, 
providing  bounds  of  the  dt^  and  dz^^  are  known. 


II- 4  Relation  Between  Time  and  Distanne  Errors 

To.-furi^r  simplify  matters,  an  approximate  relation  exists  between 

^  Md  .  “Partial  “differentiation 'of  "Eq(  4)  yields  the  two  expressions 

OZk 

(after'“’‘transposing.  -and  nrultipiyihg...by.  'Z~^ ) ; 


boc 


=  Z 


-1  h'c 


(9) 
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and 


i  ■  •“  li.  -  S.  => 


From  the  definition  of  T,  the  i  th  component  of  the  column  vector 
is  found  directly  to  he 


St 

'St, 


=  z. 


(10) 

St 

(U) 


Differentiating  the  Z  matrix  with  respect  to  yields  a  similar  square 
symmetric  matrix  with  elements  given  by 


^  (i  ^  j) 


ij 


(12) 


and  it  is  clear  that  the  (O^O)  element  is  zero  for  any  z,  .  Then 
a  is  clearly  a  column  vector  with  coD^onent  value 


'Sz 

icf 


J=0 


a . 
J 


(15) 


St 

Moreover,  is  a  column  vector  whose  i  th  component  is 
^^k 


St 


=  1 


i-1 


(ih) 


By  rearranging  and  combining  terms  of  (15)  and  (14),  the  component  value 
^  /St  SZ  — \  , 

'Sz  "  ^  t?6COin6S 

k  k 


St  Sz  — 

A  ■  “ 


J  1 


=  1  z. 


i-1 


<1 

J=0 


^k  -J 


This  simplifies  to 


-  z, 


ZI  (J) 

J=0  k  J 

(15) 


■  •>  J 


r  St  ‘  I  Sz  -  1  _  .  i-1  ^  i  / 


(16) 
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upon  making  use  df'Eqns.  (2)  and  (3)^  and'  letting  be  the  missile 


velocity  at  Since  ordinarily  represents  a  veiy  smaJi  time  interval, 

while  z;,/v  represents  approximately  the  time  interval  for  the  missile 
to  traverse  the  distance. between  the  k  th  station  and  the  origin,  the  ,E^ 
terms  may  be  neglected'.'  'Then  the  i  "th  component  becomes,  approximately,. 
-z^/'Vj^;  but  the  column  vector  with  these  components  is  just  (- 
by  Eq'(ll).  Therefore,  a  good  approximation  to  (lO)  is 


(TT) 


Making  use.. of,.  (9)  and.(l7),,  Etl( 8)  becomes: 


da 


N 


C 

k=l 


(18) 


Eq(l8)  demonstrates  that  an  "equivalent  timing  error"  can  be  constructed 
from  the  pure  time  error  and  the  pure  distance  error;  this  assumption  was 
made  in  reference- Ol"X  and  from  the  derivation  of  (18)  appears  to  be 
justified. 


An  expression  for  the  variance  of  the  elements  of  da  follows  from 
(18),  and,  assuming  corresponding  conditions,  yields  results  identical 
with  those  in  [jlj.  Letting  ff^(x)  represent  the  variance  of  x,  we  have 


(19) 


This  equation  (in  different  notation)  is  used  in  reference  to  discuss 
the  effect  of "the  spatial  distribution  of  observation  stations  on  errors 
in  the  a  elements,  and  an  equation  similar  to  (18)  is  employed  to  estlinhte 
upper  bounding- errors  for  a  five  station  range  in  reference  [2]. 

From  (19)>  the  variance  of  the  observation  errors  is  seen  to  be  the 
sum  of  the  variances  of  both  chronograph  and  distance  errors.  We  shall 
now  determine  the  variance  of  the  pure  time  error  to  be  expected  from 
digital  type  chronographs,  in  order  to  compare  it  with  the  distance  error. 
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III.  CHRONOGRAPH  ERROR 


III-l.  Model  Chronographs 

To  find  a  mathematical  expression  for  the  chronograph  error,  we 
consider  model  chronographs  with  properties  which  correspond  generally 
to' those’ of  'act'ual  chronographs  .  "We  "assume : 

(a)  All  chronographs  used  have  identical' properties; 

(■b)  The  driving  signal  from  the  oscillator  to  each  counter  is 
periodic,  with  period  e; 

(c)  The  driving  signal  is  admitted  to  the  digital  counter  by  a 
gate  circuit  controlled  by  identical  start  and  stop  signals; 
however,  the  effective  response  time  of  the  chronograph  to  the 
start  signal  is  allowed  to  differ  from  that  of  the  stop;:  signal. 

III-2.  Time  Observations 

The  oscillator  pulses  provide  the  time  scale  by  dividing  time  into 
adjacent  intervals  each  of  duration  e.  The  leading  edge  of  each  interval 
(that  part  of  the  interval  at  smallest  time)  is  numbered  consecutively. 

In  addition  to  the  e  interval,  we  consider  two  other  time  intervals,  a 
and  p,  which  depend  on  the  response  of  the  chronograph  to  start  and  stop 
signals,  respectively.  These  times  are  made  definite  by  the  following 
definitions : 

(a)  If  the  start  signaT'-begins  at  a  time  t^  within  the  interval 

.  L  -  .  given  byL-,  (.  ,  q. 

me  .-i-  .a  Z.  ^  .(m+l)e  +  a  ,  (20) 

then  the  m+1  interval  causes "the "first  count  to  be  recorded 
in  the  counter; 

(b)  If  the  gate  closing  signal  begins  at  a  time  t,-,  within  the 

JIj 

interval  given  by 

ne.  +  p  Z.  tg  Z  (n+l)e  +  p  ,  (21) 

then  the  n  th  interval  is  the  last  one  recorded -in  the- counter  . 

We  consider  all  possible  paired  values  of  time  at  which  the  start  and 
stop  signals  can  occur  within  the  above  limits  to  give  a  counter  reading 
(in  e  units)  of  n-m.  From  (20 )  and  (2l),  the  maximum  time  interval  is 
r(n+l)e  +  pj  -  jme  +  aj  ,  while  the  minimum  is  |ne  +  pi  -  jjm+l)e  +  a]. 
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The  difference  between  the  maximum  and  minimum  intervals  is  2e, 
and  the  mean  interval  is  (nvm)e  -  (a-p).  The  quantity  7  (=a-p) 
is  a  time  interval  representing  the  differing  response  of  the  chronograph 
to  start  and  stop  signals,  and  may  be  either  positive,  negative,  or 
zero.  In  an  actual  chronograph  the  magnitude  of  7  is  probably  much 
smaller  than  e;  but  we  retain  it  in  our  analysis  for  the  sake  of 
generality.  I 


The  reading  on  the  k  th  coimter  is 
true  time  interval,  ~  only  tc 

a  stopping  error  and  7  as  follov 


\  =  ^k  -  ^ 


thus  seen  to  represent  a 
within  a  starting  error. 


s; 


+  5 


2k 


(22) 


The  errors  5^^^  can  each  have  any  val\|i 


e  between  +  e/2  with  the  same 

probability;  the  distribution  for  this  ty^e  of  function  is  therefore 

2 

rectangular,  with  a  variance,  a  {&.,),  given  by 

JK: 


e^/12 


The  time  data  used  in  the  column  vector  t 
readings  converted  to  suitable  units .  Th^ 
counter  readings  in  lieu  of  the  "true"  tir.i 
is  found  that  the  arrangement  and  coupling 
bank  affects  the  size  of  the  error  in  the 


III-5.  Arrangement  of  Counters;  N  Counteas 


It  is  necessary  to  distinguish  betweq 
counters.  If  all  counters  are  driven  by 
said  to  be  "phased";  if  they  are  driven ’bjf 
of  the  same  frequency,  the  counters  are  s 


(25) 

of  Eq(i;)  are  the  counter 
error  incurred  by  using  the 
les  will  now  be  examined.  It 
of  the  counters  in  the  counter 
important  elements  of  a. 


n  two  methods  of  coupling  the 
ne  frequency  source,  they  are 
separate  unphased  sources 
did  to  be  " unphased". 


The  first  arrangement  considered  is  one  with  N  counters  to  in¬ 
dicate  time  values  at  the  N  observation  stations;  all  counters  receive 
a  simultaneous  start  signal  prior  to  the  instant  of  the  first  observation. 
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and  are  stopped  in  seq.uence  vlth  the  corresponding  .distance  observations 
at  each  station.  We  assume  the  "true"  tijne  origin  to  he  that  of  the 
first  ohservation,  and  the  time  at  the  k  th  station  is  t'  .L  where 


+  ^21.)  -  t^LL  *  "21] 


c 


The  difference  in  the  readings  of  the  first  and  the  k  th  counters  repre¬ 
sents  the  true  time  to  within  the  indicated  errors.  With  this  choice 
of  time  origin  there  appears  to  he  no  error  at  the  first  station,  although 
the  timing  error  at  any  other  station  could  he  as  large  as  2e.  However, 
we  can  take  advantage  of  the  results  of  paragraph  II-2,  and  use  a  new 
time  scale,  t,  where  t  =  t*  +  give 

\  =  <=1.  -  <=1  ^  ^  •  (25) 


Now  represents  the  "true"  time  to  within  the  error  (6^  +  ^  e 

for  every  station  including  the  first. 

If  the  counters  are  " unphased",  the  variance  of  the  timing  error 
2  2 

a  ^2k^  each  station  is  e  /6.  But  if  the  counters  are  phased, 

8^  is  constant  independent  of  k.,  Again  moving  the  "true"  time  origin, 
in  this  instance  to  we  can  replace  (25)  hy 

‘i*'  '  (26) 

2 

the  variance  of  the  time  error  is  now  €  /l2.  The  use  of  phased  counters 
has  thus  halved  the  variance  for  this  counter  arrangement. 

III-4.  N-1  Counters 

An  alternative  counter  arrangement  makes  use  of  N-l  counters  for 
N  ohservation  stations;  all  counters  are  started  simultaneously  with  the 
first  ohservation  and  stopped  in  sequence  at  the  succeeding  stations. 

The  couinters  are  numbered  from  2  to  N.  Assuming  the  origin  of  the  "true" 
time  t'.,  to  coincide  with  the  first  observation,  the  time  of  the  k  th 
ohservation  is 
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For  unptiased  counters,  -y  represents  the  "trUe®'  time  at  all  stations, 
save  the  first,  to  within  the  error  there  is  no  timing  error 

at  the  first  station.  Thus,  comparing  the  unphased  counter  arrangements, 
it  appears  that  slightly  "better  data  is  o"btained  with  the  N-1  arrangement 
than  with  the  N  arrangement'  "because  there  is  no  timing  error  at  the  first 
station  when  N-1  counters  are  used.  For  this  accirracy  to  be  obtained,  it 
is  essential  that  the  value  of  y  be  known. 

If  phased  counters  are  used  in  the  W-1  arrangement,  (=  has 

the  same  value  for  any  k^  redefining  the  time  origin  by  t  =  t*  -  6 
enables  us  to  replace  (27)!  by 

\=Cj^-7  +  62k 


and 


Then  using  -  y  for  the  time  at  each  station  except  the  first,  where 
the  value  6  is  used.  Implies  that  there  is  an  error  with  magnitude  be- 
tween  the  values  +  e/2  at  every  station;  the  error  variance  is  e  /l2. 

These  eriors  are  comparable  to  those  of  the  N  phased  counter  arrangement. 

Althoiagh  the  discussion  of  errors  implies  that  data  from  N-1 
counters  is  at  least  as  good  as  that  from  N  counters,  the  N-1  arrange¬ 
ment  has  serious  practical  disadvantages.  Aside  from  the  fact  that  the 
value  of  y  must  be  known,  we  note  "that  a  failiire  to  start  the  N-l 
counters  at  the  first  station  resiUts  In  the  loss  of  time  data  at  every 
station.  For  this  reason  alone  the  additional  counter  may  be  worth-while. 
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IV.  COMPAEiaoW  OF  TIME  AIID  DISTANCE  .EBRORS 


The  effective  ohservation  error  is  made  up  of  hoth  time  and  distance 
error  according  to  (19)*  Since  these  errors  are  of  parallel  importance, 
their  comparative  sizes  are  useful  in  deciding  on  the  chronograph  reso¬ 
lution  time  needed.  For  example,  a  hank  of  phased  1.6  megacycle  chrono- 

graphs  (e  =  .625  n  sec)  yields  a  timing  error  variance  of  €  /l2;  assuming 

-3 

the  r.m.s.  distance  error  to  he  10  feet  at  each  station,  the  separate 
contrlhutions  of  time  and  distance  error  to  the  total  error  are  the  same 
at  a  missile  velocity  of  ft. /sec.  At  lower  velocities  the  distance 

error  is  the  larger  one. 

Similarly,  to  make  the  contrlhutions  the  same  at  a  missile  velocity 
of  10,000  ft. /sec.  using  10  megacycle  phased  counters,  the  distance 
error  should  he  reduced  to  .0035  inches.  Considering  the  difficulties 
of  measuring  distances  on  the  order  of  50  feet  to  within  an  error  smaller 
than  this  makes  it  appear  that  little  advantage  is  gained  hy  using 
counters  with  smaller  resolving  times  unless  ultra-high  velocities  are 
contemplated . 

Since  the  reproducihility  of  the  experiments  depends  also  on 
physical  quantities  other  than  the  time-distance  history  (such  as  the 
fluid  medium  and  the  similarity  of  models),  errors  in  these  other  quanti¬ 
ties  should  he  considered  in  efforts  to  make  measurements  of  greater 
precision. 

G.  D.  KAHL 

F.  D.  HEHNETT 
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The  precision  of  the  least  squares  method  of  fitting  time-distance  data 
from  on  observed  trajectory  is  considered  for  the  case  vhere  tijue  Is  assumed 
to  be  a  polynomial  in  distance.  A  relation  between  the  variances  of  time  and 
distance  errors,  assumed  in  a  previous  work,-  is  derived.  With  this  relation 
and  a  model  of  digital  chronographs,  the  variance  of  the  timing  error  is  found 
for  selected  chronograph  arrangements j  in  a  favorable  arrangement,  this 
variance  is  only  one-twelfth  as  large  as  the  squared  counter  resolving  time. 

A  camparlson  of  time  and  distance  errors  is  made. 
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The  precision  of  the  least  squares  method  of  fitting  time-distance  data 
from  an  observed  trajectory  is  considered  for  the  case  where  time  is  assumed 
to  be  a  polynomial  in  distance .  A  relation  between  the  variances  of  time  and 
distance  errors,  assumed  in  a  previous  work>  Is  derived.  With  this  relation 
and  a  model  of  digital  chronographs,  the  variance  of  the  timing  error  is  found 
for  selected  chronograph  arrangements;  in  a  favorable  arrangement,  this 
variance  is  only  one-twelfth  as  large  as  the  squared  counter  resolving  time. 

A  camparison  of  time  and  distance  errors  is  made. 


The  precision  of  the  least  sqtoares  method  of  fitting  time-distance  data 
from  an  observed  trajectory  is  considered  for  the  case  where  time  is  assumed 
to  be  a  polynomial  in  distance.  A  relation  between  the  variances  of  time  and 
distance  errors,  assumed  In  a  previous  work,  Is  derived.  With  this  relation 
and  a  model  of  digital  chronographs,  the  variance  of  the  timing  error  la  found 
for  selected  chronograph  arrangements;  in  a  favorable  arrangement,  this 
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to  be  a  polynomial  in  distance.  A  relation  between  the  variances  of  time  and 
distance  errors,  assianed  in  a  previous  work,  is  derived.  With  this  relation 
and  a  model  of  digital  chronographs,  the  variance  of  the  timing  error  is  found 
for  selected  chronograph  arrangements;  In  a  favorable  arrangement,  this 
variance  is  only  one- twelfth  as  large  as  the  squared  counter  resolving  time. 
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The  precision  of  the  least  squares  method  of  fitting  time-dletance  data 
from  an  observed  trajectory  is  considered  for  the  case  where  time  is  assumed 
to  be  a  polynomial  in  distance .  A  relation  between  the  variances  of  time  and 
distance  errors,  assumed  in  a  previous  work,  is  derived.  With  this  relation 
and  a  model  of  digital  chronographs,  the  variance  of  the  timing  error  is  found 
for  selected  chronograph  arrangements;  In  a  favorable  arrangement,  this 
variance  is  only  one- twelfth  as  large  as  the  squared  counter  resolTlsg  time* 

A  ccmparlBon  of  time  and  distance  errors  Is  made. 


